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I. CONTROL SYSTEM OVERVIEW 


This thesis is intended to act as the preliminary ground 
work for the attitude control system of the Naval 
Postgraduate School satellite ORION. The background for 
ORION is detailed in Reference 1. The serbittal control 
system deals with the motion of the center of mass in an 
orbit around the earth, while the attitude control system 
deals with the angular motion of a satellite around its 
center of mass. The attitude control system must 
incorporate attitude sensors, control laws, actuators or 
thrusters, the effect of disturbing forces and the dynamics 
of the satellite. 

ORION is being designed to be launched via the extended 
Get-Away-Special (GAS ) cannister aboard the Space 
Transportation System (Space Shuttle). ORION is currently 
under development and many of the systems which will 
determine its final characteristics have not been fixed. 
For the purpose of this thesis, ORION is considered to be a 
cylindrical satellite of uniform mass with dimensions of 
0.889 meters height and 0.4826 meters diameter. rt nhasewea 
mass of 113 kilograms. Figure 1.1 shows a diagram of ORION. 

The design of a control system is highly dependent on 
the mission of the satellite. Reference 2 provides a 


section on the analysis of mission requirements and the 
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Figure 1.1 ORION Diagram 
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possible control systems which will fulfill mission 
requirements. Tables 1 and 2 show the trade-offs that need 
to be made when choosing a type of control system ([Ref. 
2:pp. IV-12,IV-12]. 

ORION is being designed to act as a payload bus for a 
variety of missions. General requirements are that the 
satellite fly in orbits from 200 kilometers to 2000 
kilometers. Attitude sensor accuracies will not be less 
than 2.0 degrees. The life of the satellite will depend on 
a mission payload versus fuel trade-off. For missions that 
do not require a high degree of pointing accuracy, spin 
stabilization provides the simplest means of attitude 
control. A spinning satellite will be less susceptible to 
nutation if spun about a maximum moment of inertia and will 
be stable if energy can be dissipated. The equations 
governing satellite motion are developed in Chapter III. 

The purpose of the control system is to maintain the 
satellite at a specified orientation in space. If the 
satellite is perturbed from this orientation by external 
factors such as gravity gradient torques, aerodynamic drag 
or nutation then the control system must sense the change 
and apply the control torques necessary to reorient the 
satellite. It must also be able to reorient the satellite 
to accommodate any change in the mission requirements. 
Figure 1.2 shows a diagram of the basic control system. The 


feedback loop and error determination process may occur 
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Figure 1.2 Attitude Control Block Diagram 
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either on board the satellite or via a satellite to earth 
station telemetry link. On-board processing requires 
preprogramming the reference angles into the processor and 
1s limited to those periods when the control system needs to 
act only as a regulator and any disturbances are small. A 
telemetry link 1s required to control the satellite during 
major orientation periods and flight maneuvers. TiteeeaalisG 
acts as a back-up for the on board processor. 

There are three distinct phases of the attitude control 
problem. They are: 1) launch, consisting of activities 
required to take the satellite into a preliminary orbit; 2) 
aeqguisition, in which the satellite's preliminary 
Orientation and maneuvers are accomplished; and 3) mission 
operations, in which the mission requirements are carried 


Siivpeet. 3S:p. 12]. 


A. LAUNCH 

Control during launch is limited to the launch vehicle 
and has little relationship to the attitude control of the 
satellite. Of primary concern is the detachment of the 
satellite from the launch vehicle in some preliminary orbit. 
The satellite will have the angular momentum of the launch 


vehicle plus some velocity imparted from the separation 


process. ORION will be carried into orbit by the Space 
Shuttle and launched from a GAS cannister. Reference 1 
details the cannister and the launch process. The launch 


Ibis 


vector is along an outward radial from Earth's center at 4 


feet per second. 


B. ATTITUDE ACQUISITION 

The acquisition phase commences after ORION is detached 
from the Shuttle. During this phase the control problem 
consists of determining the attitude of the satellite, 
spinning the satellite up to the desired spin rate, 
reorienting the satellite to align the maneuvering thruster 
with the velocity vector required to place the satellite in 
its mission orbit and maintaining the correct attitude 
GQuring the orbit transfer. Once in its final orbit, the 
control system must reorient the satellite to the desired 
mission attitude. 

NASA prohibits the firing of thrusters over 5 pounds 
sooner than 45 minutes after separation. Thrusters under 5 
pounds cannot be fired closer than 200 feet to the shuttle. 
At a launch velocity of 4 feet per second, this means the 
satellite will drift for 50 seconds before beginning spin up 
for attitude orientation. During this drift time the 
primary forces acting on the satellite are aerodynamic drag 
and gravity gradient torque. Additional torques will be 
present if the satellite contacts the cannister side during 
the launch process. The environmental torques are developed 
in Chapter III. 

The attitude of the satellite can be predicted during 


the drift period by applying these torques to the basic 
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satellite equations of motion. Some assumptions must be 
made concerning the initial conditions at the time of 
launch. These are: 

(1) The satellite does not impact the cannister during 
the launch process. If this happens the satellite 
will tumble and the attitude cannot be predicted. 
The motion would have to be observed by the Shuttle 
crew and reported. 

(2) The satellite is aligned as shown in Figure 1.3 with 
the z axis pointing towards the earth and the x axis 
in the direction of flight. The initial values of 9, 
Yand © are zero. 

(3) The aerodynamic drag force acts only against the 
satellite's center of pressure and opposite to the 
Grrectiommof flight. 

After the satellite is 200 feet from the shuttle, it is 
spun up to 60 rpm and the actual attitude determined by the 
horizon and sun sensors. A description of the horizon and 
sun sensors is given in Chapter II. Once the orientation of 
ORION is determined, control torques may be applied to 
orient it to any desired attitude. During this time system 
tests may be performed to verify the correct operation of 
the satellite. Details concerning control torques are found 
ine enapter Lil. 

The next portion of the acquisition phase is to orient 
ORION for orbital maneuvering. It is unlikely that the 
Shuttle will be able to launch ORION into its final orbit. 
For its first mission ORION will be required to extend four 
magnetometers from the body. This deployment will have two 


major effects. It will despin the satellite and it will 


increase the moments of inertia about the spin and 
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Figure 1.3 Initial Axes Orientation and Attitude Angles 
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transverse axes. Once the magnetometers are deployed, the 
satellite will need to be oriented to align the propulsion 
thruster with the velocity vector which will maneuver ORION 
mee, 1tS final orbit. Details on the final orbit are not 


yet available. 


Cc. MISSION 

The third phase of the attitude problem is the control 
of the satellite while it is performing its mission. 
Generally this requires the control system to act as a 
requmator to Maintain the correct orientation. If the 
required attitude is fixed, the control loop may be 
programmed into a micro-processor and function autonomously. 
The telemetry link would be used to monitor the processor 
and act as a back-up. Control by the earth station may be 
necessary to orient the satellite for maneuvers or to 
correct for a large disturbance. Figure 1.4 summarizes the 


Semerol proplem an algorithmic form. 


D. ASSUMPTIONS 

The purpose of this thesis is to form a model of ORION 
and simulate its motion during the three phases of the 
control problem. It will model the satellite's equations of 
motion, the environmental disturbance torques, and the 
control torques and study the interactions between these 
factors. As stated earlier, ORION is still in the early 


stages of design. Many of the systems which determine the 
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flight characteristics have not been developed. he 1s 
necessary therefore to make a number of assumptions which 
simplify the simulation problen. As the design of ORION 
matures assumptions may be removed and the model updated. 
This will help in determining design trade-offs and increase 
the accuracy of the model. 

The first major assumption is that the satellite has a 
uniformly distributed mass. This locates the center of mass 
at the origin of the body aligned axis as shown in Figure 
1.1. The moments of inertia are dependent on the deployment 
of the magnetometers, but in the launch configuration I, is 
3.239 kg-m? and I} is 8.916 kg-m*. As the mass distribution 
changes because of component location or changes in fuel 
weight, payload, etc., these moments of inertia will change. 
If the moments of inertia are located along the body aligned 
axes, the cross products of inertia are zero. The second 
major assumption is that the orbit will be circular. This 
results in a constant angular velocity 4%, around the earth. 
The radius of orbit also remains constant throughout the 
Oisbit . 

The model of the satellite will first be constructed 
assuming the satellite is a rigid body with no internal 
energy effects such as sloshing fuel tanks. It will also be 
assumed that there will be no significant changes in the 


moments of inertia because of the control thruster burns. 
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Energy effects can be taken into account later in the 
simulation process. 

The model will also assume that the spin rate of the 
satellite will remain constant. In actuality the drag force 
on the satellite will cause the spin rate to decrease 
significantly in approximately 50 days [Ref. 4:p. 34]. This 
will be ignored in the model since if a_ significant 
reduction is noticed, then the spin thrusters can easily be 
fired to spin the satellite without affecting the 
orientation. 

When a satellite is in orbit, it is subjected to a 
variety of environmental forces and moments. These 
environmental disturbances include gravity gradient torque, 
aerodynamic drag, solar radiation pressure, magnetic field 
torques and micrometeor and radiation impacts. Gravity 
gradient and aerodynamic drag effects are dependent on the 
radius of the orbit. As the radius increases these effects 
are decreased. Solar pressure and magnetic torques are 
important over a long period of time, but can be ignored 
when looking at the short term motions. The primary effect 
of these environmental disturbances iS noticed in the 
nutation and precession of the spinning satellite over time. 
This model will assume that the total torque on the 
satellite will be the sum of gravity gradient torque, 


aerodynamic drag torque and applied control torques. The 
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model will allow any additional torques to be applied at the 


Simulator's discretion. 
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fi. ATTITUDESSENSOKS 


Before control can be applied to a satellite, the actual 
orientation in space must be determined. This is done by 
the attitude sensor package. Many devices and methods’ can 
be used as attitude sensors. These include rate gyroscopes, 
star sensors, sun sensors, earth limb or horizon sensors, 
and magnetometers. This chapter will describe how 
magnetometers, sun sensors and horizon sensors operate. 
These three types of sensors are currently planned for use 


in ORION. 


A. SUN SENSORS 

Sun sensors are used to measure the angle of the sun 
relative to the spin axis of the satellite. The sensor 
consists of a pattern of photoelectric cells and a housing 
that limits the angle at which the sun's rays contact the 
cells. Figure 2.1 shows a schematic of the sun sensor [Ref. 
3:p. 163]. 

The field of view slit is oriented parallel to the x-z 
plane and perpendicular to the spin axis. The angle 6 de- 
termines the amount of sunlight on the photo cells which in 
turn produces a given amount of voltage. As the spin axis 
moves in inertial space, usually due to a torque, & changes 
and the voltage from the photo cells changes. The accuracy 


of the sun sensor is affected by the quality of solar cells 
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Figure 2.1 Sun Sensor Schematic 


27 


used, the presence of any highly reflective components of 
the satellite and the location of the sensor on the body. A 
typical accuracy is 0.01 degrees [{Ref. 3:p. 17]. The rate 
of change of 8 is determined by calculating the time rate of 
change of the output voltage. 

Sun sensor data can be used to calculate nutation and to 
calculate the yaw angle ’. The method for determining 
nutation is described later in this chapter. ¥Y is deter- 
mined by using information from the sun sensor and the 
horizon sensors. The motion of the spin axis is a combina- 
tion of “both iGeand... If 6 is known from the horizon 


sensors then for small angles, 

y2 = g2 = 62 (25m) 
and 

Y . 3a (2.2) 


B. HORIZON SENSORS 

Earth limb or horizon sensors are IR sensors with a set 
field of view that measures the difference in temperature 
between space and the earth's limb. Space is considered to 
have a temperature of 4 degrees Kelvin and the Earth has a 
temperature of 300 degrees Kelvin. As the sensor's field of 
view sweeps from space onto the earth a pulse is started. 


The pulse ends as the FOV goes from earth back onto space. 
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Two sensors are located equidistant from the center line 
of the satellite as shown in Figure 2.2. As the satellite 
spins, the two sensors create a pair of pulses as they cross 


the earth. The pulse width is determined by 


PW = R/2hé (253) 
where R is the earth chord length (Figure 2.3), h is the 
altitude of the satellite and 6 is the spin rate of the 
satellite. 

The sensors can determine d, 0, and 8. The spin rate 6 
can be computed by measuring the time between successive 
pulse starts. ¢% is determined by comparing the pulse widths 
of sensor 1 and sensor 2. When the pulse width of sensor 1 
is less than the pulse width of sensor 2, ¢$ is positive. 
When the pulse width of sensor 1 is greater than that of 
sensor 2, ¢ is negative. This relationship is shown in 
Pagure. 2.3. The value of 6 is computed by comparing the 
magnitudes of the pulse widths with those of 6 = 0 and using 
the methods in Chapter 11 of Reference 5. 6 is determined 
by comparing the pulse widths of successive pulses over 
time. 

The accuracy of the sensor is dependent on the quality 
of components, and the presence of infrared sources other 
than the earth. Masks are used to block the effects of the 
sun and avoid the saturation of the sensors. The satellite 


itself may act as a noise source, however shielding and 
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Figure 2.2 Horizon Sensor Location 
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Figure 2.3 Relationship Between Earth Chord and Roll Angle 


on 


knowing the infrared signature can reduce the effects. 
Pages 365-409 of Reference 5 give computational methods used 
to eliminate errors and determine the accuracy of the 
sensors. In general the accuracy of a horizon sensor can be 


assumed to be 0.1 degrees [Ref. 3:p. 17]. 


C. MAGNETOMETERS 

Considering the accuracy in attitude measurements 
required for ORION, the sun and horizon sensors will suffice 
for the sensor package. Since magnetometers are going to be 
used in the mission, they can be used to provide additional 
attitude inputs. These inputs can be compared with those 
obtained from the sun and horizon sensors to help eliminate 
noise and errors in meaSurements. 

For attitude sensing purposes, three mutually orthogonal 
magnetometers are required. These may either be internally 
mounted or deployed on booms [{Ref. 3:p. 250]. The 
magnetometers consist of a coil around a ferromagnetic core 
of known permeability and area. When in the presence of a 
magnetic field, a voltage is produced that is proportional 


to the field component along the solenoid axis and is given 


by 
V = -ANU(dGB;/dt) (2) 


where B; is the field component [Ref. 3:p. 181]. By 
combining the field components from the three magnetometers, 


the solenoid axis directions can be computed. If the 


32 


solenoid axes are aligned with the body axes, the attitude 
angles can easily be found. 

Errors in the magnetometers may be the result of the 
model used to predict the earth's magnetic field, the 
influence OF the satellite's electrical components, 
crosstalk between the magnetometers or a misalignment of the 
solenoid axis. At high orbital heights, the exact magnitude 
and direction of the earth's magnetic field is not well 
known and the model used to predict the magnitude and 
Girection of the flux lines may include considerable errors. 
At altitudes above 1000 km, the field strength decreases as 
the inverse cube of the distance and the magnetic field of 
the satellite may dominate. Crosstalk involves’ the 
influence of one magnetometer on the others as a result of 
eddy currents that are produced. Misalignment may occur if 
the solenoid axis is not aligned with the body axis. This 
may occur upon deployment or if there is some deformation of 
the boom during a spinning maneuver. [Ref. Sie sar 


181,182,250] 


D. NUTATION 

All sensor measurements described so far have assumed 
that the satellite was not nutating. If nutation occurs, 
the sensor data will exhibit an oscillatory behavior. The 
frequency of oscillation is the same as the body nutation 


rate wp). This motion can be used to calculate the nutation 


33 


angle directly from sensor data. The derivation and details 
are given on pages 539-548 of Reference 5. 
If the nutation is monitored by a sun sensor, then the 


maximum nutation angle iS given by 


n = 68/2R, (2.6) 


where 68 1S the range of nominal sun vectors and Re is the 
ratio of observed sun angle variations to the maximum 
nutation amplitude. R, is dependent on the satellite's 
moments of inertia and the location of the sun sensor [Ref. 
37D. go oOlr 

Knowing the nutation angle and the attitude angles, the 
Satellite can determine when to fire control thrusters to 
reduce nutation and reorient itself. This can either be 
done autonomously or via the telemetry link. Examples of 


active nutation control and orientation are given in Chapter 


Ve 
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IIIf. ATTITUDE DYNAMICS AND CONTROL 


Be ALTITUDE DYNAMICS 

Newton and Euler first described the motion of a system 
of particles in inertial and rotational frames of reference. 
These are described and detailed in a number of 
publications. This chapter will present a summary of these 
derivations as they relate to modelling ORION's motion. 

The satellite operates with respect to a reference frame 
which rotates around the earth in the orbital plane. As was 
shown in Figure 1.3, the x axis is in the direction of 
flight, the z axis points towards the earth and the y axis 
is perpendicular to the orbital plane. 

The orbital axes are related to the fixed reference 


frame of the earth by the direction cosine matrix Cc. Cisa 


linear combination of Euler rotations defined as 


where Fy is the body reference frame and /; is the inertial 


reference frame [Ref. 5:p. 23]. The three rotations are 
1 0 0 
C(o) = 0 cos ¢ sin ¢ (321) 
0 -sin ¢ cos 6 | 


ys) 


c(0) = 0 1 0 (30en 
sino 0 Cosmem 
cos ¥ Sin ¥ 0 

C(y) = -sin ¥ cos ¥ 0 (ares) 
O O 1h 


The most commonly used sequence is the 3-2-1 sequence which 


results in 


cos¥cos@ sinYcos®¢ sin¥ sind 
-sin!’ cosOsino +cos!? sir cosé 


Cc = -sin’¥cosQ cosY¥cos@¢ cos¥sino (oe) 
+sinY sindsin®g ~cosh Sir sind 
-sind -singcosd coso0cos¢ 
The satellite orientation is related to the orbital axes 
by the body fixed coordinate axes (x,y,z). The deviations 


of the body axes and the orbital axes are given by the 
angles 4 (phi), Y (psi), and O (theta). 
There are six Euler equations that govern satellite 


motion [Ref. 5:p. 32]. They are: 


p = wtp + f (3.5) 
Pp = mv (3.6) 
Re = Vv (Si7) 
h = -wtho + T (see 
h = Iw (35.9) 
e. = Sete (Zen) 
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Where Gee is defined as 


1 
0 — Wy Wy 
y* = Wo ) Wy (3%) 
Wy Wy 0 


Equations (3.5), (3.6), and (3.7) are the translation 
equations of motion and are not discussed in this thesis. 
Equations (3.8), (3.9), and (3.10) apply to the attitude 
Gomerol of ORION. 

The moment of inertia about an axis is the sum of the 
masses times the square of the distance of the mass from the 


axis. The moment of inertia matrix I is given as 


I= Iyx Ty Iyz oe 2) 


where Iy = {x2 dm, ie f y? dm, I, = f22 dm. Since ORION 
is a symmetrical body with the body axes centered at the 
center of mass, the cross-products of inertia, Iyys Ivz and 
Ioy are zero. 


Expanding Equation (3.8) into x, y, and zZ components 


yields 
hy = Why - %hz + Ty (Ses) 
hy = Wyhy - Yehy + Ty (3.14) 
hz, = Wyhy - Why + Tz (Sess) 


a, 


Expanding Equation (3.9) and substituting into (3.13), 


(3.14), and (3.15) and solving forale aelac 


Ty = I,wy + (IZ-Ly) wywz (3.16) 
Ty = Iyoy + (Iy Iz) wxwe (351) 
Te, = Igqzt (Ty-Iy) wywx (Sas ) 


T; is the sum of she external Conquesson™ the cartel eae omen 
the ith axis. These torques include environmental torques, 
control torques and any other torques which may affect the 
satellite. Torques will be discussed later in the chapter. 


Since ORION is symmetric about the x and z axes, I, = Iz, 


= I;, the transverse moment of inertia. The y axis is also 

the spin axis and Iy = Ia. The angular velocity vectorw 

decomposes into its Euler angle components: [Ref. 4:p. 30] 
Wy = 9 - Th, (3 5 
Wy = Cae. (3.20) 
Wo = Y + GW ( 3.6920) 


ORION will spin such that the spin axis is in the negative y 


direction. This means that 


Wy = Rage (32225) 


In the absence of external torques the motion equations 


reduce to 


re (Ta-It) wywz (373233) 


Ttuy =e (3.24) 
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It Ws = (Wem a) oy ox (3729) 


Equation (3.24) implies that Wy is a constant. Differen- 


tiating Equations (3.23) and (3.25) and substituting for 0, 


and Wy yields 


(I_-I,) (1-1 ) 


al 
tes i = 5.2.20 
We + iD WD Wye 0 ( ) 
© 
Cie eee) 
fee 2 tee tf, (Sci) 
Z Z Var’ 
he 


For stability, the coefficients of wy, and w, must be 
positive. The y axis is the maximum moment of inertia when 
I, > It and is the minimum moment of inertia when Izy > Ig. 
If there is no energy dissipation, the satellite is stable 
if spun about either the maximum or minimum moment of 
inertia. If there are external torques or energy 
dissipation means, the satellite is stable only if spun 
about the maximum moment of inertia [Ref. 4:p. 112]. 

Since there are external torques present, these must be 
taken into account when calculating the attitude of the 
Baeeoleite. Substituting for »™ in Equations (3.16), (3.17), 


emer ( 3.13), 


Te ($-UgY¥) + (Te-Za) (97 0) (Yt ty?) 


Ty 


oe * e 2 ee ° 
Ted-Ipewol + (Ie -Ta) ( Wolt Uo 9-OF — 6 09) 


ityuthes (ite it) (a Gomsou en (mete ised) aeemtent! (3228) 


Sie, 


te (3520) 


Tz = It(¥+u9d)+(La-Tt) (wo 6) (6 -wo¥) 
ee a e 2? ee ° 
= Tt¥tIeuopt (La lt) (wo ¢7wol = Ootwo OY) 
= It¥+(Ia It) (wood we) Yt (La-Tt) (wo @) +1tWe) 4 


= Tey (Ip-Ta) (same 0) Ua aC rea ea, rire ed (3.30) 


Solving for the second derivatives of the Euler angles 


yields 

6 = (Ty-(1e-Ia) (wS-uge) 6-1 ee mee 
Yo = (B=(I,-Taq) (uS-ug 6) YC Cee Ta (Ug=O) —T pug OL meen 
9 = -Ty/Ia (3.33) 


If wo and 9 are assumed to be constant the equations can be 


written 
6 = Ty - AdO- BY (3.34) 
vy = T, - Ay+ Bg (3.35) 
0 = -Ty/la (3. 36) 
where: 
A= (1 = Pe7Lie oC) ain 
ea bak Ia/It) (wo70) - Wo 


The block diagram form of the equations is given in Figure 
Beolee As can be seen from Figure 3.1, ¢ and ¥ are cross- 


coupled through their rotation rates. 
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System Block Diagram 


Benes sok 
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The total angular momentum vector H can be expressed as 


Ho = Iguy + Tpog Cee) 
2. ae Z ; ; 
where wt = Wy + we and lies in the x-z plane. Wy rotates 


relative to the satellite with the angular velocity . where 
A = (Ia/Ie = 1)wy (3.38) 
The body nutation frequency is defined as 
Wy = Ta/lIte wy (273.9) 


The motion of the satellite consists of the satellite 
rotating about its spin axis and the spin axis rotating 
about H with the angular velocity w, [Ref. 4:p. 115]. The 


nutation angle n is defined as 


IT, w 
oe ‘4 (3.40) 


Figure 3.2 shows the relationship between +, wy, wy and H. 
The inherent stability of the satellite can be found by 
analyzing the transfer functions of the satellite. Setting 


T, and T, equal to zero and dividing by I; yields 


0 (3.41) 


) +(1-I4/Ie) (wa 0) Wop Ay eee ery 


Y(t) (1,5 - 6) UY een ste (W—-0) wo) 4 0 (2425 


Noting that -\ = (1-92/714) (@.5-¢°) and) taking (Chee tap rae 


transforms of the equations yields 
6 (s)s* - Wod(S) - Atwo)s¥(s) = 0 (235 
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43 


¥(s)s* - Awe ¥(s) + (Atug) So(s) = 0 (3.44) 


Solving for Y¥(s) and substituting 


ee 
5 (Atm) S 
0(S).(S —AQ))) = ——— a Ss) — 0 G45) 
(s -Aw_) 
O 
oS) (S2-Aw)? + (At+W9)282O(s) = 0 (3.46) 
6(S) (S4+(Q +Wo) 2-lWo) 8% + A2W5) = 0 (3.47) 
Partitioning 
BD 
AW 
s O (3.48) 
O'S) tac’! ae ee 
s©(s°+((Atw,) “-Aw,)) 
The open loop transfer function is 
yD 
a (3.49) 


2 ee 2 

s (s + ((A+W,) ~hw,) 
The transfer function has two poles at the origin and two 
poles on the imaginary axis at +j V(A+u,)2-Aw, and 


-3/(, +49) 2 Awo- Therefore, the stability of the satellite 


depends on spin rate, inertia ratio and orbit altitude. 
Table 3.1 shows the root locations for different o and Wy « 
Because there is no cubed term in the transfer function, 
the roots will always be on the imaginary axis. This 
implies that the satellite will tend to reach a new steady- 
State value for ¢and Y if perturbed, rather than return to 


an initial state. 
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B. ENVIRONMENTAL TORQUES 


Before the motion of ORION can be accurately modelled, 


the possible external disturbances must be studied. 


The 


major disturbances experienced by a satellite in low earth 


orbit are aerodynamic drag and gravity gradient torques. 


Aerodynamic drag 


satellite through the upper atmosphere. 


as 


where: 


D 


is a result of the motion of the 


It can be expressed 


= OV Cy A (3.50) 


density of air at orbit altitude 


velocity of satellite along orbit 


drag coefficient 


surface area over which the drag acts. 
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The drag on the body can be represented by a Single force 
acting on the body's center of pressure. If this center of 
pressure is located away from the center of mass, a torque 
is developed. 


The velocity of the satellite is 


V = oR (3.51) 
2 
WR ya (3.52) 
v2 = ,./R = ee (3.53) 
where: 
Ug = gravitational constant 
R = distance of satellite from earth's center. 
The density of air decreases as the altitude increases. At 


a 200 km orbit the air density is approximately 7.0 x10719 
kg/m>. The surface area is the cross-sectional area of the 
Satellite; for ORION, A = 0.429 m*. Drag coefficients vary 
with shape, material smoothness, etc. A typical value for 


this type of satellite is Cp = 2. The drag can be expressed 


as 
14 
np = Lixo 6 (3.54) 
R 
The center of pressure for ORION is assumed to be 
located 0.076 m along the positive y axis. Since the drag 


acts opposite to the direction of flight, the total drag 


torque on ORION can be expressed as 
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i = > eecos 7 (3.55) 


where L is 0.076 m. 
The second major disturbance is due to the gravity 
gradient that exists as the distance from the earth 


increases. The gravity gradient torque is defined as 


3u 
G = == he xR) (r-R) dm (3.56) 
R 
O 
where: 
R = position vector of the satellite with respect 


to the earth; and 


r = coordinate of an elemental mass from the 
Cenreryonr Mass. ihRef. 4:p.. 131) 


Using the 3-2-1 Euler rotation between the body axes and the 


orbit fixed axis, the gravity gradient torques are 


Gy, = BOT (test -) sine cos ? cos~0 } (32a57 ) 
Gy = 6 (372355)) 
Gp = 306 (Ip-Iq) sin6 cos 0 sin $} (Seo) 


Thus gravity gradient torques depend on orbital altitude, 
moments of inertia and satellite orientation. 

The equations describing aerodynamic drag and gravity 
gradient torques can now be used in simulating their effect 
on satellite motion. The analysis of this effect is 


contained in Chapter IV. 
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C. COMPUTER MODELS 

The computer model was written in the Dynamic Simulation 
Language (DSL). The program listing is contained in the 
Appendix. The program was intended to be as general as 
possible to allow it to apply to a variety of parameter 
changes. Of special interest are I+, Ig, orbital altitude, 
air density, and spin rate. 


The three Euler equations are integrated once to obtain 


od, VY and ©. These are integrated again to obtain 6, ¥, and 
ee The initial values depend on the _ situation being 
analyzed. Once the attitude angles have been calculated, 


they can be used to calculate “;, “,, A =and the nutation 
angle. 


Also included is the equation governing magnetometer 


deployment. It is implemented in the form of a ramp 
FUN POTN. The starting time is the time when the 
Magnetometer is to be deployed. The ending time is 


calculated by dividing the desired boom length by 0.0169. 
If no boom length change is desired the ramp function can be 
commented out of the program and the variable Ly, set to the 
desired boom length. An Ly of zero indicates that the 
magnetometers are in the storage position. 

Changes in spin rates can be accomplished in two ways. 
One is to change the initial condition of © in a PARAM 
statement. This is useful when evaluating the effect of 


configuration changes at different spin rates. The second 
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method simulates the use of the spin-up thrusters. Ty can 
be a step, ramp, or pulse input, or any combination that 
approximates the actual thruster firing profile. 

In addition to the variables describing moments of 
inertia, spin rate, and thruster torque, many other 
components of the model are left as variable quantities in 
the program. These can be changed by PARAM statements in 
multiple runs to provide the widest range of possible 


parameter studies. 
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IV. PARAMETER STUDIES AND ACQUISITION SIMULATION 


A. SIMULATION MODEL 

The design of any modern system involves the 
investigation and analysis of how the system responds to 
inputs with a given set of system parameters. These 
parameter studies are used in design and cost trade-offs 
which will eventually lead to the optimum system design. 
One way to run parameter studies is to use scale models and 
testing chambers. These are primarily used in the final 
design stages. In preliminary design where there are a lot 
of variables computer simulations are best performed. 

The computer model described in Chapter III was used to 
analyze some parameter changes that may be associated with 
ORION. It is expected that this model will be used to 
determine the optimum configuration in terms of moments and 
magnetometer deployment. This chapter analyzes the changes 
in inertia ratio and spin rates as a result of a typical 
mission. It also analyzes the environmental effects of 


gravity gradient and aerodynamic drag. 


B. PARAMETER STUDIES 

The first study was to analyze how a change in the 
moments of inertia would affect the spin rate and stability. 
The current mission requires the deployment of four 0.907 kg 


(2 1b) magnetometers on lightweight booms. The booms were 
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assumed to deploy at a constant rate of 0.0169 m/sec to a 
maximum of 3.048 m (10 ft). The mass of the booms was 
considered to be negligible and did not contribute to the 
inertia of the satellite. The magnetometers are located in 
the x-z plane of the satellite as shown in Figure 4.1. 

The effect of deploying the magnetometers is to increase 
Born I, and TI. The new moments of inertia can be 
calculated by finding the moment of inertia for the 
satellite without the magnetometers and then adding the 


moment of inertia for each magnetometer. Thus 


I, = O0.5mor? + 4my(lph + r)2 (4.1) 
oe —MeOmeSema(sr-th-) +e2m,,(1p+r) - (ez) 
where: 
Ms = mass of satellite minus the mass of 4 
magnetometers 

Mn, = mass of one magnetometer 

1g = radius of the satellite 

h = height of the satellite 

lh = distance of the magnetometer from the outside 


Opeunesmcatel lite. 


Substituting in the values for ORION, the moments of inertia 


can be found at any boom length by 


BeeOcmemceo2ay ly, +0), 242) (4.3) 


Ta 


3.8m 1.91(1,, + 0.242) - (4.4) 


It 


Sl 





MAGNETOMETER 


Figure 4.1 Magnetometer Deployment Diagram 
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Figure 4.2 shows this relationship between the boom length 
and the moments of inertia. I, becomes larger than I+ at a 
boom length of 1.53 m. If energy could be dissipated, then 
with I, larger than I+ the satellite would become stable 
about the spin axis [Ref. 4:p. 117]. 

Deployment of magnetometers also has an effect on the 
spin rate and could be used to change it without using the 
Spin-up thrusters. Using the principle of the conservation 


of angular momentum, if wt is assumed to be zero then 
H = Tay (4.5) 


and as I, increases Wy must decrease. Figure 4.3 shows how 
the spin rate changes with boom length. 

As noted earlier, the magnetometer deployment can be 
used to change the spin rate of the satellite. If the sat- 
ellite is spinning at 60 rpm (Wy, = -6.283 rad/sec), a boom 
length of 1.097 m will slow the spin to 20 rpm (Oe =2-094 
rad/sec). A 1.87 m deployment will slow the spin to 10 rpm 
(wy = —1.047 rad/sec). Adjusting the boom length to reach a 
specific spin rate is applicable if the mission does not re- 
quire the magnetometers a specific distance from the satel- 
lite. It is more probable that the boom must be a specific 
length and the satellite must spin at a specific rate. 

There are two ways to achieve this configuration. The 
first method is to deploy the magnetometer to the desired 


length, then adjust the spin as necessary with the 
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Change in Moments of Inertia with Boom Length 
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Figure 4.3 Change in Spin Rate with Magnetometer Deployment 
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thrusters. A second method is to despin the satellite with 
the thrusters then deploy the magnetometers and readjust the 
spin with the thrusters. If the change in spin and boom 
length is small, it may be possible to despin the satellite 
to a rate such that the subsequent deployment of the booms 
will despin the satellite to the required rate. 

Since the firing interval of the thrusters is 
proportional to the amount of fuel used, it can be used to 
compare the two methods. To illustrate this, assume the 
Satellite is initially spinning at 60 rpm with I, = 3.239 
kg-m*. A boom length of 2 m and a spin rate of 20 rpm is 
required. At 2 m, I, is 21.447 kg-m?. Using the first 
method, after deployment the spin rate is 


ta 
ie a il 
yy a2 


(3° 232 ie occ 3) 
ZAG 


-0.9476 rad/sec (4.6) 


Assuming the total torque of the thrusters is FxR = 0.2168 


n-m, the firing interval dt is 


T7 (2.094 —W) 
F XR 


) 
dt = Ye 


(2 Aaa 00 4 Oa) 
0.2168 


= [isa 3sesec (4.7) 
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Using the second method, the firing interval required to 


despin the satellite to 20 rpm is 


(Sez eon Greco] 2. 093) 


dt foeolle oS) 


62.6 sec (4.8) 


Deploying the magnetometer to 2 meters causes the spin rate 
to decrease to 0.3152 rad/sec. The firing interval required 
to readjust the spin is 176.1 sec. hae Wee caleeet lieing 
interval for the second method is 238.7 sec, more than twice 
that of the first method. 

This indicates that the first method is more fuel 
efficient in deploying the magnetometers and achieving a 
specific spin rate. A disadvantage of both these methods is 
that if the booms are too flexible, adjusting the spin may 
cause the magnetometers to lag the boom-satellite attachment 
points and cause too great a stress in the booms. 

ie ~ Stability 

Whether © is changed by deploying magnetometers or 
by changing the internal mass distribution of the satellite, 
it will have an effect on the stability of the satellite. 
Figure 4.4 shows how a change ino affects 4 and &,. Since 
the roots of the open loop transfer function depend oni, 
they also change with o as was shown in Table 3.1. Figure 
4.5 shows the responses of the satellite to a 1.0 n-m pulse 
of 0.01 sec as o changes. This was done at three different 


spin rates of 60 rpm, 20 rpm, and 10 rpm. The advantage of 
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Figure 4.4 The Effects of Inertia Ratio on 
Nutation Frequencies 


aS 


8.3 9.9 
Time (SEC) 


t 6.7 6.6 6.9 1.8 


9.6 6.1 6.2 4.3 6.4 4.5 6, 
Sire. (SEC). - ——_—_——_— sates 


et eee 4.2 9.3 0.0 4.8 .e 
Tim€ (S€C} 


Figure 4.5 Satellite Response to a Short Duration Torque 
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increasing 0 1s to make the satellite less susceptible to 
nutation. Figure 4.6 shows that the higher the spin rate 
and the larger the inertia ratio, the smaller the nutation 


will be for a given input torque. 


C. ENVIRONMENTAL EFFECTS 

The aerodynamic drag and gravity gradient torques were 
modelled using the equations in Chapter III. These 
equations were incorporated into the simulation to see what 
effect they would have on the satellite. 

Since the drag force is dependent on the density of the 
atmosphere, it decreases as the orbit height increases. 
This is shown in Table 4.1. At an altitude of 200 km the 
drag torque is 1.383 x1073 n-m and decreases to 1.762 x107? 


at an altitude of 1000 km. 


TABLE 4.1 


TYPICAL VALUES FOR AERODYNAMIC DRAG 


ALT (Km) Density (kg/m?) Drag (n) Tp (n-m) 
100 6 x1077 15.84 1.204 
200 7 x10719 .0182 1383 <0 
400 1 x1071t1 2.523 x 1074 1.917 >eom. 
1000 1x 10719 2.318 x 10m 1276221 


Gravity gradient torque is also dependent on orbital 


altitude since it depends on wo. Unlike drag, however, 


60 


2.0 


1.6 


1.2 





NUTATION ANGLE (DEC) 
0.8 


0.4 


0.0 


0.30 0.46 0.80 0.76 0.90 1.06 1.20 1.35 1.69 
INERTIA RATIO (SIGMA) 


Figure 4.6 Nutation Angles for Different Spins and Inertias 
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gravity gradient torque is not a constant but depends on the 
moments of inertia and the orientation of the satellite. At 
an altitude of 200 km the maximum torque is 3.99 X1079 n-m. 
This decreases to 3.72 10715 n-m at 1000 km. 

An orbital altitude of 200 km was used to determine the 
effect environmental forces have on ORION. Inertia ratios 
of 0.364 and 1.488 and spin rates of 10, 20 and 60 rpm were 
used to give the widest ranges of results. These results 
are shown in Figures 4.7, 4.8 and 4.9. At 10 rpm the 
maximum nutation is 2.053 1073 radians for a o of 0.364 and 
5.33 107° radmans for a Gd of Meat ocemen cece rpm the nutation 
angles are 5.173x 1074 radians and 1.17 x107> radians and 


for 60 rpm they are 6.0*x 10°72 radians and 1.33) xtoee 


radians. 
AS can be seen in the figures, the nutation is 
oscillatory. This is because of the oscillatory nature of 


the environmental torques, principally the gravity gradient. 
When the satellite nutates to its maximum angle, the torque 
effect is reversed causing the nutation to decrease. When 
the nutation is zero, the torque effect iS again reversed 
and the process repeats itself. The frequency of the effect 
is related to the spin rate of the satellite. This is seen 
in the 0 term of the torque equations. 

If the accuracy requirement of the attitude control 
system is small enough, environmental effects must be taken 


into account in the simulations. They will be ignored in 
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Figure 4.7 Nutation Due to Environmental Torques; 10 Rpm 
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Figure 4.8 Nutation Due to Environmental Torques; 20 Rpm 
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Figure 4.9 Nutation Due to Environmental Torques; 60 Rpm 
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further simulations described in this thesis since their 
effect is less than the 2 degree accuracy requirement of 
ORION. This will simplify the model without a large error 


in accuracy. 


D. ACQUISITION PHASE PROFILE 
After analyzing the effects of oO changes and 
environmental torques, the next step was to simulate the 
motion of the satellite from launch to spin-up. Drag and 
gravity gradients effects were included in this simulation. 
The flight profile consisted of launching the satellite 
with no spin. All initial values of ¢, ¥ and 0 were assumed 


to be zero. At 50 seconds a step input of Ty = 0.2168 n-m 


waS applied. This caused the value of Wy to increase. At 
144 seconds the input was removed. The value of Wy was 
-6.29 rad/sec (60 rpm). Figure 4.10 shows howo and ¥ 


change during this portion of the mission. 

During the drift period environmental torques cause the 
satellite to rotate 0.013 degrees in roll and 11.04 degrees 
in yaw. As the spin increases the o and ¥ begin an 
oscillatory motion. The amplitude of the motion depends on 
the spin rate. As the spin rate increases, the satellite 
becomes more resistive to torque effects. When the spin 
rate becomes a constant, the amplitude of oscillation due to 
environmental torques becomes a constant although the scale 
in Figure 4.10 is too large to show this. At 144 seconds, ¢ 


is 3.89 degrees and ¥Y¥ is 12.54 degrees. The nutation 
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360 
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caused by environmental torques is 0.0015 degrees, much less 
than the acceptable nutation of 2 degrees. After 400 
seconds with no controlling inputs, ? is 7.4 degrees, ¥ is 
10.8 degrees and the nutation is 0.0015 degrees. 

From this point, the magnetometers could be deployed or 
control torques could be applied to reorient the satellite. 
Chapter V discusses the methods of controlling nutation and 


reorienting the satellite. 
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V. MISSION PHASE SIMULATION 


Pee OAL bLITE CONTROL 

The process of maintaining the correct orientation of a 
satellite consists of eliminating the precession and 
nutation effects caused by external or internal 
disturbances. The simplest case is when the spin axis and 
angular momentum vector are aligned and are perpendicular to 
the orbital plane. 

This presents a two part control problem. The first 
part is to reduce nutation by forcing u, equal to zero. The 
second part is maintaining the angle between the satellite's 
angular momentum vector H and the orbital angular momentum 


vector Hj. The relationship between H and Ho, is 
heer 
O 
Cos a a 
lal HI 
The angle & can be obtained by precessing the satellite with 
control torques and then eliminating any resulting ;. 
W, can be forced to zero by applying a control torque of 
magnitude equal to Iywy in the direction of Ip-¥e. % can be 


calculated from measurements One, to) > and %. The 


magnitude of torque required is determined by 
fe oe edt == TeWt 


where: 
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F = the impulsive force of the thruster 


r = distance of the thruster from the Cencer 
of mass 
dt = time over which the thruster is fired. 


An ideal thruster would provide an impulsive force at the 


correct time in the spin cycle and dt would be Zero. 


The 


firing profile of an actual thruster can be approximated by 


the trapezoid model shown in Figure 5.1 [{Ref. 3:p. 21}. 


is the time the thruster is commanded on, ty is the 
thruster reaches maximum force, ty is the time the 
is commanded off and ts 1s the time the thruster 


zero. The total time dt is from t, to tes. 


at which the effective torque is applied. 


ty 


time the 
thruster 


force is 


te is the time 


tp and ts can be determined from experimental or 


manufacturer's data. Knowing t,, ty and ts, t, and ty 


be determined by 





Tan Ato = bDfa (5. 
Gye 
a acos At +bsin At (5. 
ia Cc ®: 
where: 
cos Ats - cos Aty cos Ate - cos At, 
a = UC) ee oe 
A(t, t,) A(t. ty) 
. sin At, - sin At) - sin Ate - sin At, 
= ——Y\(t,-=) 0.6 )° 
A(t, tq) A(t. ty 


70 


the 


can 


1) 


2) 


ross 


Figure 5.1 Typical Thruster Profile (Trapezoid Model) 


fick 


Figure 5.2 shows the relationships between the firing 
profile and timing. 

An even simpler approximation is the pulse function. 
This would occur if ty = t, and ty = ts. Knowing that dt = 


Ca-to, tos (17 27Gb sae Therefore 


to t. — (lone (5.3) 


th to GaGiyiiat (545) 


Once the magnitude of control torque and the firing arc 
are known, the problem is when during the spin cycle to fire 
the thruster to eliminate %;,. The acceleration of any point 


(x,y,z) on the satellite is 


ay = =xWS + wnYA sin jt C5)..5) 
ay = (ug-A)A(Z cos At + x Sin \t) (5:26) 
az = ohne + wnyAcos At (537) 


where A iS proportional to the nutation angle [Ref. 4:p. 
AEG ae If an accelerometer is placed in the x-z plane, the 
axial acceleration reaches a positive peak when %; points 
towards the accelerometer. Placing a thruster 90 degrees 
from the accelerometer produces a torque opposite to wy. If 
the accelerometer is located on the Z axis, 

ay = (Wg7\)AZ cos At C5..op 


and the peak acceleration is reached when day/dt = 0. 
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Figure 5.2 Relationship Between Spin Rate and 
lea ielalig(s; Ubsbuvar 


ES 


Since 
day/dt = -(Wg-A)AZA sin At (5.9) 
then 


O when t = TA 


day/dt 


The value of acceleration at the firing point is 


Using the value of t = tf as t,, the on and off times for 
the thruster can be calculated. If the spin rate is too 


fast, dt may encompass several revolutions of the satellite. 
To provide a torque in a given direction, however, the 
firing must be Jimited™ to Jan arc centered abcuraeee 
Therefore, the total firing time required (dt) must be 
divided into segments occurring at integer multiples of tg. 
For example, assume the nutation angle is 1.0 deg 
(0.0174 rad), Izy is 8.91 kg-m?, wy is -6.283 rad/sec, and H 


is 20.35 kg-m?/sec*. Therefore 


TeWe H sinn 


0.355 kg-m?/sec¢ (5.11) 


Assuming F is 0.448 n and R is 0.242 m 





dt = => 


3.28 sec . G5i,.1 2) 
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At 60 rpm the satellite makes one revolution per second. If 


the firing arc is 90 degrees, the firing time per arc is 


i = 0.25 sec (oO - 13)) 


Ng 


Therefore, 13 firings will be required to eliminate the 
nutation. If the spin rate was 10 rpm, only 2 firings would 
be required. Each firing must occur at aymax to avoid 
adding energy to the system. 

This method of active nutation control was simulated 
using the computer model. Nutation was started by inputting 
a 5 n-m pulse of 0.01 seconds. After 0.03 seconds the 
control torque could be applied. The minimum required 
nutation was set at 0.001 radians to avoid computational 
errors. The firing arc was limited to 45 degrees. 

To apply the correct control torque, the program 
calculated the value of ay which would cause _ the 
acceleration peak to be in the center of the firing 
interval. When the calculated value of ay was greater than 
this value, Tz, was set equal to 2Fr/dt where dt is the 
firing arc time period. As long as the nutation was greater 
than 0.001, T, was applied at the peak positive acceleration 
point. 

Mis simasation was run for 3 different spin rates of 
10, 20 and 60 rpm at an inertia ratio of 0.363. It was also 
run for two different inertia ratios and a constant spin 


rate. These results are shown in 3 series of 3 figures. 
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Each series shows the nutation angle and applied control 
torque, the ¢and Y angles and the firing diagrams for a 
different spin rate. 

In Figure 5.3, the spin rate was 10 rpm. MThe initial 
nutation was 1.451072 radians. It was reduced to 0.001 
radians at 4.4 seconds. Tz, waS applied to reduce the 
nutation to 0.001 radians. From Figure 5.4, 6 and Y were 
reduced to steady-state values of 1.67 <1072 radians anda 
-1.66 x 1074 radians. The values of 6 and ¥ oscillated about 
these values because the nutation was not completely 
eliminated. Figure 5.5 shows the firing diagram for this 
case. Ty, was 8.33 x 1072 n-m over a firing time of .6 sec. 


In Figure 5.6, the spin rate was changed to 20 rpm which 


resulted in an initial nutation of 7.28x 1073. The first of 
3 control firings occurred “aeeZz. 2 ceconds. Nutation was 
reduced after 12 seconds. In Figure 5.7, the values of9 
and y are 7.91 x1072 radians and -2.1 x1074 radians. From 


the firing diagram of Figure 5.8, Ty, was 4.1-x 1072 n-m over 
a pulse of .36 sec. Three pulses were needed for a total 
firing time of 1.08 seconds. 

In Figure 5.9, a spin rate of 60 rpm was used. The 
Init ral nutation was PAs On radians, the first 
application of control occurred at .7 seconds and the 
nutation was reduced after 27.5 seconds. g and ¥ from 
Figure 5.10 were 2.27x1073 and -2.0 x107> radians. From 


the firing diagram of Figure 5.11, 18 applications of Ty, 
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Figure 5.3 
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32 


40 


were required. T, was 1.3 x107% n-m with a pulse width of 
.23 seconds. A total firing time of 4.2 seconds was 
required. 

Figure 5.12 shows the difference in control for the same 
spin rate of 20 rpm but for two different inertia ratios. 
When ©O = 0.363 it required 3 applications to reduce the 
nutation. When waS increased to 1.12 the initial amount 
of nutation was less and only one application of control 
torque was required. 

As a result of these simulations, it can be seen that 
the slower spin rates require less firing times and number 
of pulses. This is because the firing arc time is larger. 
A low number of pulses also reduces the probability of a 
timing error which may cause nutation instead of eliminate 
lee In every case shown, ¢? reached a steady state value 
while Y approached zero. This is a result of using T, for 
the input and control torques. 

Fewer pulses resulted in shorter firing times. and 
quicker responses. This equates to a fuel savings. The 
model can also be run with different size thrusters by 
changing the value of F in the program. Larger firing arcs 
would also decrease the response times especially at higher 


spin rates. 


B. SATELLITE REORIENTATION 
The control system may be used to reorient the satellite 


as well as reduce the nutation. Reorientation can 
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MY 


20 Rpm 


20 


be accomplished either as part of the nutation reduction or 
as a separate control sequence. It was shown in Figures 
5.4, 5.7 and 5.10 that although the nutation was reduced, $ 
and ¥ were not returned to their initial values. 


Reorientation during nutation control is a two-step 


process. The first step involves reducing the nutation 
angle by half its original amount. This will cause a new 
nutation cone to be established. When the spin axis moves 


180° along the new cone, control torques are applied to 
reduce the nutation to the minimum level. 

Because of the differences between nutation frequency, 
spin rate and}, the thrusters may not be aligned with the 
wt vector when the satellite reaches the 180° pornt on tae 
nutation cone. The time required for the spin axis to 


travel 180° is 
At = /wn (5.14) 
The satellite will rotate about its axis through an angle 
ol Se ee (5. DoF 


in this time period. The relative rotation of the 4, vector 


is 
B= At (5.16) 
over the same time period. The two angles ®% and 8 are equal 


only Ifo —==2207 
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Figure 5.13 shows the relationship between Wy, Yy, and dh. 


Since 
= (o-l)wy 
= (Ia/It - 1) wy (S575 
and 
Wy = Ta/Tey (open?) 


The two angles a and 8 can be calculated as 


an ole (5.19) 
ee 7/6 (5.20) 


Since over the required time period the thrust vector may 
not be aligned with the ™; vector, the application of a 
control torque would add to the nutation. 

To overcome this difficulty, four sets of thrusters can 
be used. If the thrusters are located on the x and Z axes, 
a combination of thrusters can be fired to make the 
resulting torque collinear with the wy, vector. 

This is accomplished in the program by saving the time 
at which the nutation is halved as FTIME. When the program 
time reaches FTIME plus At, 8 is calculated and the correct 
combination of T, and T, used to reduce the nutation to an 
acceptable level and restore ¢and ¥ to their initial 


values. 
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Figure S23 





Relationship Between Wy, A, anduy 
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The second method is to reorient the satellite assuming 
there is no nutation. The degree of reorientation was given 


earlier in the chapter as 


Knowing the angle a and H,, the H required to make a zero 
can be calculated. This gives the required w+. A control 
torque equal to Itw; is applied to start the satellite 
nutating. When At is equal to T/w, a control torque equal 
and opposite to I;yup is applied to reduce the nutation. 
This is done ina manner similar to that of reorienting the 


Sacelinteswiath nutation. 
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VI. CONCLUSIONS 


This thesis derived the equations which describe the 
rotational motion of a satellite. It also derived the 
equations describing the effects of gravity gradient torque 
and aerodynamic drag. These equations were then used to 
construct a computer Simulation model. The model was used 
to analyze how changes in moments of inertia and spin rates 
affect the motion of the satellite. 

It was shown that extending magnetometers on booms would 
cause I, and I; to increase. This changed the inertia ratio 
o. The spin rate of the satellite could then be changed by 
changing the boom length of the magnetometer. As o 
increased, the stability of the satellite was changed. 

A high value of o made the satellite resistive to nutation 
by external forces. 

The effects of gravity gradient and aerodynamic drag 
were analyzed. The amount of nutation caused by environment 
depended on the spin rate and o. The maximum amount of 
nutation occurs at a low spin rate and lowo. The nutation 
is on the order of 2.053 x107~* radians for a o of 0.363 and 
a 10 rpm spin rate. The effects of the environment need 
only to be taken into account in the motion analysis if the 
accuracy required by the control system is small. The 


accuracy for ORION is to be within 2 degrees so gravity 
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gradient and drag were only applied during the example 
mieicghte profile. 

The control system that is used is an active nutation 
control. An example of how the control system would be used 
was given. Assuming a firing arc of 459, 4.4 seconds were 
required to reduce the nutation at a spin rate of 10 rpm, 12 
seconds for a spin rate of 20 rpm and 27.5 seconds for a 
spin rate of 60 rpm. The procedure for reorienting the 
satellite to any desired attitude was also described. 

Although the computer model was written to simulate the 
satellite as closely as possible, many areas need more 
detailed modelling. The satellite was assumed to be a 
cylinder of uniform mass. The actual mass distribution 
needs to be determined and the moments of inertia 
calculated. This would change the inertia ratio. and 
possibly cause some cross-products of inertia to develop. 

Another factor that was not included in this model is 
the effect of liquids slosh in the fuel tanks. This could 
cause internal torques that adversely affect the motion. 
These tanks could also be designed to dissipate energy. If 
this could be done, then when © was greater than one, the 
satellite could be made stable. 

The firing profile of the control thruster was assumed 
in the model to be a pulse function. Entering the profile 


as a trapezoid would increase the accuracy of the model. 
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When real thruster data becomes available it could then be 
entered. 

As has been stated before, many of the variables in the 
program can easily be changed in parameter runs. This 
allows a great deal of flexibility in analyzing design 


changes as they occur. 
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APPENDIX 
ATTITUDE ANGLE CONTROL PROGRAM 


ieee ATTITUDE ANGLE CONIROL PROGER&? 
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Peo ROGhere es nov LDES TRE SIMULATION MODEL FOR THE CONTROL OF THE 
OO eo eee eUokoe tHe BOQUATIONS DERIVED IN CHAPTER 3 TO MODEL 
SATeLLiIT& MOTION AND ENVIRONMENTAL TORQUES. IT CALCULATES THE ANGULAR 
Pee OUP NGl Ee Seee tae AlTETUDE ANGLES.” 271 CALCULATES THE MAGNITUDE OF 
ieee Peeler nbd OROUES BASED ON ORISTAL PARAMETERS. IIT ALSO 
COVPURES sl oe NUTATION >) ANGULAR MOMENTUM VECTOR AND THE CONTROL TORQUES 
Maeno oan ey Oh tase aClIVE NUZATION CONTROL SYSTEM. CHANGES [O THE 
PROGREI CANS Seeriebs VIA PARAM STATEMENTS GR DIRECTLY IN TRE LISTING. 


VERITABLE DESCRIPTION 


L =U totetieemenel Ins CENTER OF PRESSURE TO THE CENTER OF MASS 
FE =—(rmrUsoler ones OF ONE CONTROE THRUSTER 

R = RADIUS Of THE SATELLITE 

EE aol oy Shoe MOMENT OF INERTIA 

IA = HOMENP PGS INERTIA ALONG THE SPIN AXTS 

ALT ="ORB It Geebtliube OF THe SATELLITE 

RO = RADIUS OF Toe EARTA 

ROE = eNoti yes AlR Aljeeee ORBITAL ALTITUDE 

WO = ANGULAR VELOCITY OF THE SATELLITE WRI EARTH 

D = DRAG FORCE 

WA = ANGULAR VELOCITY ABOUT THE & AXIS WRI EARTH 

WZ = ANGULAR VELOCITY ABOUT THE Z AXIS WRT EARTH 

WY = ANGULAR VELOCITY ABOUT THE Y AXIS WRT EARTH 
Ealb = SATELLITE ROTATION ABOUT THE X AXIS 

Polo = SeteLbite ROTASICN ABOUT THE Z AXIS 

aD = SATELLIT=2 ROTATION ABOUT THE Y AXIS 

seul = ROLL ANGLE 

jasysl = YAW ANGLE 

TH = PITCH ANGLE 

TDA ="DRAG TOnOvE iN fhe a DIRECTION 

tDZ = UhAG, (CReUe ty tre 2 DIRECTION 

AB, = TOTAL DRAG TORQUE 

GX = GRavity GRADIENT IN THE xX DIRECTION 

GZ = GRAVITY GRADIENT IN THE Z DIRECTION 

GG = TOTAL GRAVITY GRADIENT TORQUE 

as = COl@acr.s LOMOUES IN TRE A DIRECTION 

TZ = CONTROL TORQUE IN THE 2 DIRECTION 

ite = CONTROL TORQUE IN THE Y DIRECTION 

WT = PRS VSRoe ANGUEAR VELOCITY 

HT = TRANSVERSe& COMPONENT OF THE ANGULAR MOMENTUM VECTOR 
HS = AXIAL COMPONENT OF THE ANGULAR MOMENTUM VECTOR 
i = MAGNITUDE Of THE ANGULAR MOMENTUM VECTOR 

NU = NUTATION ANGLE 

SIGMA = INERTIA RATIO 

Been Rebar ive VELOCITY OF WI WRT THE SATELLITE ROTATION 
WIT = NUTATION FREQUENCY WRI INERTIAL SPACE 

AY = AXIAL ACCELERATION OF TRE SATELLITE DURING NUTATION 
FARC = FIRING aRC OF THE CONTROL THRUSTERS 

ae SmiiPeeOr FEAR POSITIVE AAIAL ACCELERATION 

RON Seno een tONeA Lethe TRHAUSTER TURN CN TIE 
Aiea — LiMe OF THE THRUSTER FIRING 


CANCeL = TIME FOR THE SPIN AXIS TO TRAVEL PI RADIANS 


CONST IT=8.911, IA=3.239, R=0.242, F=0.448, NUREF=0.001,L=0.076 


Hz 


OD RASEA 


Pie Dp FLAG 
PeeeGeR 


* 


eee ReGmrnt olAlEMENTS: ORBIT HEIGHT, DENSITY, ORBIT VELOCITY 


Per 


BS 


ALT = 200.00 

= 7.0D-10 
= (6378.0+ALT)*1.0D03 
WO = SORT (3 .987D14/RO**3) 
1.7D14*ROE/RO 








plAGu— ao 
INCON WKO=0.0,WZO=0.0,WYO=-2.093, PHIO=0.0, PSIO=0.0, THO=0.0 
WXD = (TX+TD+GK- (TT-1h) 
WZD = (TZ+TDZ+GZ-(IA-IT)*WYAWX) /IT 
WYD = TY/IA 
PHID = WX + PSI*WO 
PSID = WZ-PHIAWO 
THD = WY+WO 
WX = INTGRL(WXO,WXD 
WZ = INTGRL(WZO,WZD 
WY = INTGRL(WYO,WYD 
PHI = By ae 
PSI = INTGRL(PSIO,PSID 
TH = INTGRL(THO, THD) 
DYNAMIC 
* ENVIRONMENTAL TORQUE CALCULATION; USE WHEN DISTRUBANCE TORQUES ARE 


* REQUIRED 
TDX = -D*L*SIN(TH)*COS(PST) 
TDZ_= D*L*COS (TH)*COS (PST) 
(TDK**2+TDZ**2) **O, 5 
3.0% (WO**2 S(TE-TA)TSIN Bet a 
3. 0*(WO*x*2)*(IT-IA)*SIN(TH) *COS (TH) *SIN(PHI ) 
(GX**2+GZ**2) *KO.5 
KRKKAKAEK 


= USE See a STATEMENTS WHEN NO DISTURBANCE TORQUES ARE REQUIRED 
T = 0. 


ee 


+ tH A 
MOQ 
QINDXS 
| Ss 


‘ae 

Eo 

N 
| 


k 
* STEP FUNCTION USED TO SET AN INITIAL NUTATION VALUE 
TZ = S.0*STEP(0.01)-5.0*SDeEGmoZ) 


* CALCULATION OF ANGULAR VELOCITIES, MOMENTUM VECTORS AND FREQUENCIES 
WE = (WKAR2 + WZA*Z)%*AQ.S 
HT = ITAWT 
HS = IA*WY 

H = (HT*A2 + HS**2)**O.5 

NU = ARSIN(HT/H) 

SUGMAl =" tay in 

LAMDA = (SIGMA-1.0)*WY 

WN = SIGMAXWY 


* CALCULATION OF THE FIRING TIME AND FIRING ARC 
AY= (WY-LAMDA) * R * COS(LANDA*TIME) 
FARC = 0.392699/ABS (WY) 
TC = 3.14159/LAMDA 
AON = (WY-LAMDA)*R*COS (LAMDA* (TC-FARC) ) 


* THIS STATEMENT ALLOWS THE NUTATION TO BE PRESENT BEFORE CONTROL OCCURS 
IF (TIME.LE.0.2)Gom moans 


CONTROL CALCULATION AND FIRING. IF FLAG IS 0 THE CONTROL SYSTEM WILL 
REDUCE THE INITIAL NUTATION BY HALF THEN FIRE TO REDUCE THE ATTITUDE 
ANGLES TO ZERO. IF THE FLAG IS 1 THE CONTROL SYSTEM WILL REDUCE THE 
NUTATION TO NUREF ONLY. 
IF ((NU-NUREF).LE.0.00001) THEN 
TX =0.0 
TZ = 0.00000 
ELSE 
IF (FLAG.£0.0) THEN 


> a a a a 
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DN =_(NU- NUREF)/2.0 + NUREF 
FLAG - 
END IF 
TF (AY.GT.AON.AND.FLAG.EQ.1)THEN 
FTIME = TIME 
TZ = 2.0*F*R*FARC 


TA = 0.0 
Beos 

TZ=0.0 

TA =0.0 
END) =s 


108 (HU. fo Eire DNU ) THEN 
CANCEL = FTIME + ABS(3.14159/WN) 
FLAG = 2 
J) NCB el BS. 
* CALCULATION OF THE CORRECT TRRUSTER COMBINATION REQUIRED TO REDUCE 
* THE AITITUDE 2NGLES TO A MINIMUM. 
i> (TIME. GE CANCEL. AND. FLAG.BO.2) THEN 
DT = CANCEL + 2*FARC 
CONT = 4*F*R*FARC 
Derg =seOlll~ Siler (CANCEL) + CONT*STEP(DT) 
ee (SIGHA. m0) tae 
a ((SIGMA- 1,0) ev 4159/SIGMA) -6.283 
S 


BETA = (SIGMA-1.0)*3.14159/SIGMA 


END IF 
IF (BETA.GT.0.0.OR.BETA.LE.-1.5708) THEN 
DELTA = ABS(BETA) 
TZ = -TCL*COS (DELTA 
TX = -TCL*SIN(DELTA 
ELSEIF (BETA.GT.-1.5708.OR.BETA.LE.-3.14159) THEN 
DELTA = 3.14159- ABS(BETA) 
TZ = TCL*COS (DELTA) 
TX = -TCL*SIN(DELTA) 
ELSEIF (BETA.GT.-3.14159.0R.BETA.LE.-4.7123) THEN 
DELTA = ABS(BETA) - 3.14159 
B2e= TCL*~COS (DELTA 
ie = Tesi ( DELTA 
ELSEIF (BETA.GT.-4.7123.OR.BETA.LE.-6.283) THEN 
DELTA = 6.283 - ABS(BETA) 
107 -TCL*COS (DELTA) 
ieee TCrasaN (DELTA) 
END IF 
END IF 
15 END IF 
CONTRL FINTIM=20.0, DELT=0.01 
SAVE .01, NU, PHI, PSI, WX, WZ, WY, SIGMA, LAMDA,WN,AY,AON, 
He uceeot Pubes PSD, THD, TH, TX, TY, WT 
GRAPH (G1,DE=TEK618) TIME(UN=SEC) PHI(UN=RAD ,LO=-3.0E-03,SC=1.5E-03). 
une =RAD ,LO=-3.0E-03,SC=1.5E-03) 
GRAPH (G3,DE=TEK618) TIME NU(UN=RADS) 
GRAPH (G4,DE=TEK618) TIME TX TZ(LO=-22E-03,SC=.0165) 
LABEL (G1,DE=TEK618) ROLL AND YAW ANGLES; WY = 20RPM 
TABEL (G2,DE=TEX618) YAW ANGLE; SIGMA=0.3635 & 1.479; WY=-6.283 
TABEL (G3,G9,DE=TEK618) NUTATION ANGLE; WY = 20 RPM 
LABEL (G4,DE=TEK618) FIRING DIAGRAM WY = 20 RPM 
END 
STOP 
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